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In 1994, Wegner [l[ introduced the method of con- 
tinuous unitary transformations (CUT) to many-body 
physics. Independently, a similar method was developed 
by Glazek and Wilson denoted as similarity renormaliza- 
tion scheme 0, 9 to be used in QED and QCD. In this 
non-perturbativc approach, a Hamiltonian is mapped to 
an effective, renormalized model by a series of infinites- 
imal unitary transformations specified by the flow equa- 
tion. This effective Hamiltonian shows (block-)diagonal 
structure and can be used to calculate ground state and 
low-energy properties as well as observables d, Q. A 
perturbative variant (P-CUT) allows for systematic high- 
order perturbation expansions Additionally, the 
CUT method provides a generic framework to derive ef- 
fective models for many-particle systems that incorpo- 
rate the relevant physics. The mapping of the Hubbard 
model to a generalized t-J model illustrates this point 
exemplarily |9h13|. 

Techniques based on CUT have been applied success- 
fully in various fields of many-particle physics as electron- 
phonon-coupling [3, [t| , spin-chains [l|| [I?} and ladders 
[l8l [l9j . anyonic excitations [2(il - [22j or non-equilibrium 
problems [J] to mention only a few. For an overview, sec 
Rcfs.[13and[2l 

The basic concept of CUT is to transform the prob- 
lem of (block-)diagonalising a Hamiltonian into to the 
solution of a differential equation for the Hamiltonian 



d,H(£) = [ v (£),H(£)} 



(1) 



commonly known as flow equation |42j |. The initial 
Hamiltonian -ff(O) is linked continuously to a unitarily 
equivalent (block-)diagonal or otherwise simpler Hamil- 
tonian H(oo) in the limit of infinite £. 

A drawback of the CUT method is the creation of more 
complex interaction in the Hamiltonian during the flow. 
To keep the set of equations finite, it is necessary to de- 
fine a truncation scheme and thereby keep only a finite 
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number of multi-particle interactions. These truncation 
errors result in unknown errors for the calculated phys- 
ical quantities. A common approach is to start with a 
strict truncation scheme, i.e., a scheme comprises only 
a limited number of terms in the Hamiltonian, and to 
justify this scheme a posteriori by comparing it to cal- 
culations using less strict truncation schemes [25]. This 
yields a pragmatically well- working method to control the 
quality of the calculated quantities. However, since this 
approach only compares different approximations, the re- 
lation to the exact result or rigorous error bounds are out 
of reach. 

In our work, we analyze the effects of truncation to 
the flow and show how the truncation can be treated in 
a mathematical framework. Thereby, we aim to quantify 
the truncation error a priori. This approach will lead 
to rigorous bounds for the ground state energy and the 
highest energy eigen-value for the Hamiltonian. 

The paper is structured as follows: After this introduc- 
tion, we give an overview of the CUT method. In Sect. 
ITU we analyze the truncation and define the truncation 
error. The formalism developed is applied to the double 
hard-core boson model in Sect. |HT] to provide a trans- 
parent, simple application. Subsequently, we discuss in 
Sect. HVl the modifications needed for the calculation of 
truncation errors in extended systems. For further il- 
lustration, we show results for the dimerized Heisenberg 
chain. Finally, a summary is given. 



I. THE CUT METHOD 
A. Homogeneous flow equation 

Probably, unitary transformations are one of the most 
widely used techniques in studies on Hamiltonians. They 
render a description of a Hamiltonian possible in a more 
appropriate basis in which the physical properties can be 
studied more easily. Most desirably, every Hamiltonian 
can be diagonalized by a certain unitary transformation. 
Unfortunately, this transformation is usually unknown. 

The basic idea of the CUT-method is not to search for 
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such a transformation in one step, but to bring the Hamil- 
tonian successively closer to a simpler shape by a series 
of infinitesimal transformation. Therefore, a continuous 
flow parameter £ is introduced that parametrizes the con- 
tinuous unitary transformation U{£). The Hamiltonian is 
considered to become a function H(£) = U{£)H^W{£) 
of this parameter. In this way, the initial (bare) Hamil- 
tonian is linked continuously by a unitary trans- 
formation to the rcnormalizcd Hamiltonian showing the 
intended structure = H(oo) in the limit of infinite 
£. By derivation with respect to £, one obtains the flow 
equation 



8 t H(£) =^U\£)H{£) 
= W),H{£)]. 



H{£)U{£) 



dW(£) 



Of 



(2a) 
(2b) 



The antihermitian generator rj of the transformation 
reads 



v (£) 



dU{£) 
d£ 



U\£) = -rf{£). 



(3) 



Equation @ is linear differential equation for the Hamil- 
tonian. We emphasize that also all intermediate Hamil- 
tonians H{£) conserve the full information of the system 
because they are only written in a different basis. 

Since the basis has changed during the flow, observ- 
ables may not be calculated directly using their bare op- 
erator O^ 1 but have also to be transformed by a similar 
flow equation 



d e O(£) = [r,(£),0(£)} 



(4) 



The transformation of observables was used first by 
Kehrein and Mielke to determine correlation functions 
for dissipative bosonic systems [1, Q . 



B. Generator schemes 

Up to here, the problem of diagonalization has only 
been recast in the form of determining an appropriate 
generator r](£). The key ingredient of the CUT- method 
is to choose the generator as manifestly antihermitian 
operator depending on the flowing Hamiltonian. We de- 
note the superoperator f\ : H(£) — > r)(£) = fj[H(£)) as 
generator scheme to distinguish between the mapping fj 
and the function rj(£). In this way, the flow equation 
for the Hamiltonian @ becomes non-linear, while the 
transformation of observables (jl]) stays linear. The gen- 
erator scheme has to be designed in a way that the flow 
equation has attractive fixed points where the Hamilto- 
nian has the desired structure. In this manner, (block- 
)diagonality can be obtained by merely integrating the 
flow equation 0, 0, [M H3 ■ 

For the first generator scheme introduced by Wegner 
[H, the Hamiltonian H{£) = H d (£) + H nd (£) has to be 



decomposed into a diagonal H d and a non-diagonal part 
H nd . The generator is defined as a commutator 

r,(£) = ffw[H(£)} = [H(£),H nd (£)} = [H d (£), H nd (£)] (5) 

of the diagonal and non-diagonal-part of the Hamil- 
tonian. One directly realizes that a vanishing non- 
diagonality yields a fixed point of the flow. The proof 
of convergence for unapproximated systems was given by 
Wegner |l| for finite matrices and extended to infinite 
systems by Dusuel and Uhrig [2?} ■ The generator de- 
couples eigen-subspaccs of different energy eigcn-values, 
but it is not able to treat degeneracies. In his original 
work concerning the n-orbital model [l[ , Wegner noticed 
divergences. He could avoid them via taking only terms 
violating the number of quasiparticles into account in 
the definition of H nd aiming at block-diagonality instead 
of diagonality. In this manner, the complexity of the a 
problem can still be reduced significantly because differ- 
ent quasiparticlc spaces can be studied separately. 

To overcome the problem of residual off-diagonality 
due to degeneracies, Mielke [26| introduced a generator 
scheme on the matrix level based on a sign function of 
index differences rjij = sign(i — j)hij that always yields a 
diagonal Hamiltonian. Independently, Knetter and Uhrig 
[H,0] developed a similar scheme which concentrates more 
generally on a quasiparticlc picture. In their approach, 
the Hamiltonian H(£) = J^ij Hj(£) is decomposed into 
different blocks of terms with respect to the number 
i of quasiparticles created and the number j of quasi- 
particles annihilated by the term. In this notation, the 
generator scheme acts as 



(0) 



In the limit of infinite £, the Hamiltonian converges to a 
block-diagonal, quasiparticle conservin g st ructure if the 
spectrum is bounded from below 0, l26ll27l] . Blocks with 
i 7^ j decay exponentially with rising £. During the flow, 
the quasiparticle spaces are ordered ascending to their 
energy eigen- value prj l28l] . This implies that the vacuum 
state, i.e., the state with j = 0, is mapped to the ground 
state of the Hamiltonian if it is not degenerated. 

A special feature of this generator scheme is that it 
strictly conserves the block-band structure of the Hamil- 
tonian du ring the flow. A similar generator was used by 
Stein [29|, l3f| in a case where the sign function was not 
necessary. In contrast to Wcgncr's generator scheme, the 
right-hand side of the flow equation for the Hamiltonian 
([2]) is only quadratic in the Hamiltonian's coefficients in- 
stead of cubic. 

A recent development in the field of generator schemes 
is the ground state generator [l9[ 



rj &s [H{£)] =H i Q {£)-H° j {£). 



(7) 



The definition resembles the particle conserving scheme, 
but it is designed to decouple only the zero quasiparticle 
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subspace of a system, i.e., the ground state if not degen- 
erated. It was introduced by Fischer, Duffe, and Uhrig to 
describe quasiparticles decays, since the picture of con- 
served renormalized quasiparticles becomes very cumber- 
some. 

Compared to the particle conserving scheme, ?7g S ex- 
hibits an enhanced numerical stability and saves com- 
putational ressources due to its fast convergence. As a 
drawback, it does not conserve the block-band structure 
as rj pc does. In this work, we will make use of both gen- 
erator schemes, rj pc and 77 gs . 



C. Truncation scheme 

One way to solve the flow equation ([2} is to parame- 
terize the Hamiltonian in second quantization by an ad- 
equate operator basis. In this way, the differential equa- 
tions for the coefficients of the operators are found by 
computing the commutator on the right hand side in @ 
and re-expressing the result again in the operator basis 
chosen. The evaluation of the commutator generates new 
many-body interaction processes that are not present in 
the initial Hamiltonian. They have to be incorporated in 
the Hamiltonian which results in an iterative calculation 
of the commutator. This proliferation of terms generi- 
cally yields an infinite number of differential equations 
which is intractible in practical applications. Four differ- 
ent strategies have evolved to obtain a closed the set of 
differential equations: 

On restricting to finite systems, the flow equation can 
be solved exactly both on the level of second quantiza- 
tion or on the level of matrix elements. As a drawback, 
this approach is subjected to the same limitation as other 
finite-size methods. Furthermore, the differential equa- 
tion system may still be large enough to require further 
approximations for practical applications. 

In some special cases, it is possible to obtain a closed 
system of equations by identifying a small expansion pa- 
rameter. In his original work Wegner investigated 
the n-orbital model. He was able to close the system of 
differential equations in the limit of infinite n. However, 
these advantageous cases do not show up in every system 
or are simply out of interest. In general, approximations 
have to be applied to the system to overcome the problem 
of proliferation of terms. 

In the method of perturbative continuous unitary trans- 
formations (P-CUT) introduced by Knetter and Uhrig 
[6|-l8| , the non-diagonality H n ^ is considered as small per- 
turbation to the diagonal part of the bare Hamilto- 
nian. In this way, the flow equation can be expanded and 
therefore used to apply perturbation theory up to very 
high orders. 

The non-perturbative approach, which we will use 
in this work, is dubbed self-similar continuous unitary 
transformations (S-CUT). Here a truncation scheme is 
defined that incorporates all terms considered to be im- 
portant to the problem while other terms are neglected. 



In view of the numerous degrees of freedom, the choice 
of an adequate truncation scheme is a non-trivial task. 
It has to respect the system's physical properties. Since 
the structure of the Hamiltonian does not change during 
the flow, it is named self-similar. As a rule of thumb, 
a calculation is considered to be the more reliable the 
more terms are included in the truncation scheme. Due 
to truncation, the modified flow equation reads 

d £ H(£) = T[ V (£),H(£)]. (8) 

Here we introduced the superoperator T which denotes 
the application of the truncation scheme. 

A natural description for many-body problems is pro- 
vided by second quantization. Since most of the low- 
energy physics can be preserved by low numbers of suit- 
able quasiparticles, it is useful to define a truncation 
scheme neglecting all terms that create or annihilate 
more than a given number of quasiparticles. We empha- 
sise that this truncation expressed second quantization 
does not imply any restriction of the Hilbert space which 
is to be considered a major advantage. 

We stress that the action of a Hamiltonian on a state 
containing n quasiparticles can be split into a sum of ir- 
reducible terms affecting at most m < n quasiparticles 
each (see Fig. [T]). Thus, the truncation of high-particle 
irreducible processes does not imply the complete ne- 
glect of matrix elements between states of high quasi- 
particlc numbers, but rather an extrapolation based on 
lower quasiparticlc irreducible processes. 

In extended systems, often additional truncation crite- 
ria have to be applied in order to close the set of differ- 
ential equations. An obvious choice for gapped systems 
with finite correlation lengt is to truncate according to 
the real-space range of the physical process generated by 
the term under study. We make use of this real-space 
truncation in Sect. IIVI 

Truncations cause quantitative errors in the calculated 
physical quantities and may possibly lead to divergences 
of the truncated flow, even though convergence has been 
proven for the untruncatcd flow. In their analysis of the 
flow equations for the Anderson model [25| , Kehrein and 
Miclkc used a strict truncation scheme that includes only 
contributions of types that are already present in the bare 
Hamiltonian. 

In a second step, they included a new contribution to 
the scheme and analyzed the new equations in order to 
assess its relevance. In their classification, terms that 
affect other matrix elements only by quantitative devia- 
tions and vanish in the limit I — > oo are denoted as irrel- 
evant, or as marginal if they converge to a finite value. 
Only terms that are able to change the behavior of the 
flow equations qualitatively, for instance that cause diver- 
gences if they are not treated properly, were considered 
to be relevant. This approach is suited to ensure the cor- 
rect qualitative behavior of the effective model derived 
from the CUT. But it does not provide a quantitative 
measure of the truncation errors. 
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Figure 1: Decomposition of the diagonal element for the three-particle state into irreducible interaction processes. By 

truncation of three-particle processes, only the three-particle irreducible term a^aLa^aaa^a^ is neglected. 



D. Symmetries 

For practical computations, exploiting symmetries of 
the Hamiltonian is generally very useful. To study Hamil- 
tonians of systems with an infinite number of sites, the 
use of translation symmetry is inevitable. Let us suppose 
that we have chosen a basis of operators to express the 
Hamiltonian. Each term in the Hamiltonian is given by 
one of these basis operators multiplied by a prefactor, its 
coefficient. 

Although the Hamiltonian is symmetric under a group 
of specific symmetry transformations, the symmetry of 
each individual basis operator may be lower. In this 
case, the coefficients of several basis operators fulfil lin- 
ear conditions to ensure that their combination is invari- 
ant under the symmetry group. By parametrizing the 
Hamiltonian in terms of linearly independent symmetric 
combinations of basis operators, the number of coeffi- 
cients to be tracked is significantly reduced saving both 
computation time and memory. Technically this can be 
done by selecting one operator as a unique representative 
from which the complete symmetric combination can be 
obtained by taking the sum ^ G over a specific subgroup 
of the symmetry group of the Hamiltonian. 

We emphasize that in general one has to clearly dis- 
tinguish between the symmetry group of the Hamilto- 
nian and the supcropcrator Y]q. Even if the Hamilto- 
nian displays a continuous symmetry such as the SU(n) 
spin rotation symmetry, the number of constraints that 
can be derived for the coefficients of the Hamiltonian's 
terms is limited. Excluding translation symmetries, this 
means that only a finite number of terms are represented 
by one representative operator. Therefore the number of 
symmetry operations to build the symmetric combina- 
tion from a single representative is often finite although 
the exploited symmetry is continuous. In summary, we 
take the superoperation ^ G as a technical tool to benefit 
from the Hamiltonian's underlying full symmetry group. 

Moreover, the precise meaning of J2g depends on the 
representative under study. As an example, a represen- 
tative that shares the whole symmetry of the Hamilto- 



nian, e.g., unity, is already identical to its correspond- 
ing symmetric combination. To this operator, acts 
as identity. The other limit is a representative operator 
that does not share any of the Hamiltonian's symmetries. 
The superoperator ^ G applied to such a representative 
has to generate the fully symmetric combination of basis 
operators. 

Since the sums J^g occur on both sides of Eq. @, the 
modified flow equation can be reduced to a representative 
expression requiring a sum for only one argument of the 
commutator. Correction factors have to be introduced 
since different representative basis operators appearing 
in the commutator give rise to different numbers of basis 
operators in the fully symmetric combination which they 
represent. Details on the implementation of symmetries 
in S-CUT can be found in Ref. G3 @ 



II. MATHEMATICAL ANALYSIS OF 
TRUNCATION ERRORS 

A. Effects of truncation 

On truncating the flow equation, all information about 
the truncated terms is lost. It is common practice to ne- 
glect terms which are considered to be unimportant and 
to justify this a posteriori. But even if these terms are not 
subject of the intended analysis of the effective Hamilto- 
nian, their omission leads to quantitative deviations for 
the coefficients of all terms in the Hamiltonian because 
they are linked by the differential equations. We stress 
that the commutator of the generator and a truncated 
term may result in terms that comply with the trunca- 
tion scheme, i.e., that we want to compute quantitatively. 
Thus the loss of information cannot be limited to certain 
terms only. Gencrically, truncation introduces errors in 
all coefficients of the Hamiltonian. 

Because of truncation, the transformation of a Hamil- 
tonian H(£) described by the truncated flow equation 
([8]) does not need to be unitary anymore. Therefore, the 
spectrum of H(£) will be distorted during the flow. Phys- 
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ical quantities calculated based on the effective Hamilto- 
nian arc affected by finite inaccuracies. In the following 
sections III BI E. we present a formalism to bound these 
errors rigorously. 

An additional physical consequence of truncation can 
be derived for real-space truncation schemes which ne- 
glect interactions beyond a certain range d > d max (see 
Sect. II C]) . Due to the formulation in second quantization, 
the S-CUT method is capable to handle infinite systems. 
Nevertheless, the truncation by range affects correlations 
on larger length scales. We observed that the coefficients 
of representatives in a truncated infinite system and a 
truncated periodic system with a certain size I > Lfi n 
share the same set of differential equations, if the alge- 
bra is local. Any possibility to observe that the system 
size is actually finite is masked by the truncation scheme 
if the system size is at least Lfi n = 3e? max + 1 . 

Therefore also the intensive physical properties of an 
infinite system determined using S-CUT with truncation 
range d max are identical to those of a finite system with 
a certain size I > Lfi n and periodic boundary conditions. 
The quantity Lfi n can be understood as an effective size 
introduced by the real-space truncation scheme. The 
mathematical derivation including a numerical verifica- 
tion is given in Appendix [A] 



B. Splitting the flow equation 

To isolate the truncation error of S-CUT, we start from 
the full flow equation 

d e H(£) = [r,'(£),H(£)] (9) 

for the unitarily transformed Hamiltonian H(£) with 
H{£ = 0) = H^°\ The generator rf appears instead 
of r\ because in practice the generator is determined from 
the truncated Hamiltonian H' and not from H , see Eq. 

(K2D- 

We decompose H into two parts: the solution H 1 from 
a truncated flow equation and the difference H" = H—H 1 
between the truncated and the non-truncated calcula- 
tion. Next, the flow equation © can be split into the 
system 

d t H'{£)= T[r,'{i),H'{£)] (10a) 

d t H"{£) = (1 - f)[r,'{t), H'{1)\ + W(£),H"(£)} (10b) 

of differential equations for H' and H" . As initial condi- 
tions, we choose 

JT'(O) = ff (b) (11a) 

H"(0) = 0. (lib) 

Obviously, the sum of the equations (fTU)) with the initial 
condition (fTTj) reproduce the flow equation (J9]) with its 
initial condition. 



Up to now, the generator rf is not specified. Using the 
generator scheme rj, we define the generator 

r,'(£) = ff[H'(£)] (12) 

as a function of the truncated Hamiltonian. Note that 
this choice does not violate the unitarity of the transfor- 
mation because the generator rf continues to be mani- 
festly antihermitian. 

Equation (jlOaj) provides a closed set of differential 
equations for the coefficients of the truncated Hamilto- 
nian which can be treated by numerical integration. This 
leads to an effective Hamiltonian H'(oo) — H'^ with a 
structure determined by the chosen generator scheme. 

In contrast, the full Hamiltonian H is transformed by 
a true unitary transformation, but does not need to have 
any special structure in the limit of infinite £, since it is 
transformed like an observable by rf . However, it is to 
be expected that it is close to H' if truncation errors are 
small. 

The difference H" stores the complete 'non-unitarity' 
of the transformation of H'(£). Mathematically, Eq. 
(|10b[) describes a transformation of H" via a flow equa- 
tion with an additional inhomogeneity 

K (£) = (l-T)[ V '(£),H'(£)} (13) 

depending on H'(£). This natural emergence of an inho- 
mogeneous flow equation is quite remarkable and has not 
been observed before to our knowledge. We emphasise 
that the number of equations defining k(£) remains finite 
if H' and thus rf are restricted by the truncation scheme 
to a finite number of terms. Hence the computation of 
k(£) is indeed feasible. Of course, this is not true for H". 

C. Inhomogeneous flow equation 

To solve the inhomogeneous flow equation (|10b|) . we 
use the ansatz 

H"{£) = U{£)A{l)tf{£) (14) 

with A(0) = H"(0) = 0. The unitary transformation 
U(£) is linked to the generator rj'(£) of the transformation 
by Eq. ©. The formal solution for U(£) using the £- 
ordering operator C reads 

U{£) =£exp (j v'{l')dl^j . (15) 

Using variation of parameters 

diH"{£) = [r}'(£),H"(£)) + U(£)d £ A(£)U\£) (16a) 
±[r,'(£),H"(£)] + K(£), (16b) 
leads to the equation 

e 

A{£) = A(0) + J U\£')k{£')U{£')M' . (17) 
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Therefore, the formal solution of the inhomogeneous flow 
equation (jlObp is given by 



i')di' ] 



+ J W(£')K(£')U(£')d£' J U*(£). 

(18) 

This expression (fT8")) has a very direct interpretation: All 
contributions of the inhomogeneity up to the given value 
of I are re-transformed to I = and summed. This sum 
is evaluated after a unitary transformation to the consid- 
ered flow parameter. 

D. Truncation error 



The formal solution (JT8J) enables the calculation of the 
distortion of unitarity by the truncated calculation. All 
effects of the truncation are stored in H"(£). Certainly, 
a direct calculation is neither practical nor desirable, be- 
cause it is equivalent to an untruncated calculation. But 
for the derviation of a bound of the truncation error only 
a small part of the information is essential. To assess 
the quality of the truncation, we are interested in the 
norm of H"(£). In particular, we want to focus on norms 
that are unitarily invariant, i.e., invariant under unitary 
transformations [441] . 

We apply the norm to Eq. ([T5]) . Due to its unitary 
invariance, we obtain 



\H" 



U\£')n{£')U{£')d£' 



(19) 



In addition, we used H"(0) = from Eq. (fTTb|) . To 
avoid the complicated integration of an operator- valued 
function, we apply the triangle inequality to the Riemann 
integral arriving at the upper bound 



H"(£)\\< \\uHe')K(£')U(£')\\d£' (20a) 



Mi')\\de' =-. A H (e), (20b) 



where again the unitary invariance was used. We define 
the derived quantity Ajj{£) as truncation error of the 
transformation. By construction, it is an upper bound 
for the distance between H' and H measured by the se- 
lected norm. We emphasize that Ah{£) is a scalar func- 
tion which depends only on the norm of the truncated 
terms as defined in Eq. (|13|) . It starts at zero and in- 
creases monotonically with the flow parameter £. 

Because of the finite and constant number of terms 
complying with the truncation scheme, the number of 
contributions to the inhomogeneity k(£) stays also con- 
stant during the flow. The coefficients of the terms in 



k(£) can be calculated as functions of the coefficients of 
H'{£) already known by numerical integration. This is 
the key simplification compared to the practically im- 
possible direct calculation of H"{£) or H(£). 

In the above analysis, the necessary ingredients are the 
flow equation and the truncation scheme. Hence all con- 
siderations can also be carried over to the transformation 
of observables. The truncation error of an observable O 
can be estimated analogously by 



(l-T)[ V '(l'),0'(£')) d£'>\\0"(£)\\ 

(21) 



where O is decomposed in O' and O" in analogy to (TIT]) . 
The generator is defined by the numerically accessible 
truncated Hamiltonian H' . 



E. Rigorous bounds for observables 

The truncation error Aq is a property of the entire 
transformation of O' that quantifies the loss of accuracy 
by truncation. It is desirable to have rigorous bounds for 
the accuracy of physical quantities calculated by trun- 
cated CUTs. Indeed, it is possible to obtain such bounds 
by calculating the truncation error defined by the spec- 
tral norm 



max EV (At A). 



(22) 



For hermitian operators, the spectral norm is identical to 
the maximum absolute eigen-value. 

We denote the lowest eigen-value of O' by fi' ,„ and the 

o *i mm 

associate eigen-state by \ip'). For the untruncated observ- 
able O we use f2 m i n and \ip). Since O is transformed by 
a unitary transformation, f2 m i n does not change during 
the flow whereas Q' min is changed due to truncation, for 
illustration see Fig. [SJ 

The spectral norm of O" fullfills 

\\0"(£)\\>{0"(£)) ne) (23a) 
= (0(l)) rw -(G'(£))iP'(£). (23b) 



By condition, f2 m i n is a lower bound for ( 0{£) )i/>'(£\- It 
follows 



10" 



> a, 



n' min (£) =: An min (£). 



Analogously, one obtains the inequality 

\\0"(£)\\>-{0"(£)) m 

= -(0(£) ) m + (0'(£) ) m . 
n min >n' . (t) 



(24) 



(25a) 
(25b) 
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Since ft' min (£) is a lower bound for ( 0'(£) 
tain 



4>W 



\\o"(£)\\ > -n mhl + n' min (£) = -Ao min (f). 

In summary, the truncation error 

ko(t)>\\0"(£)\\>\An min (£)\ 



we ob- 



(26) 



(27) 



is an upper bound of the deviation of the minimal eigen- 
value of the effective operator due to the truncation. 
Analogously, one can prove that Aq defines an up- 
per bound for the deviation of the maximal eigen-value 

A very useful result ensues by considering the special 
case of the truncation error of the Hamiltonian itself be- 
cause the ground state energy can directly be read off 
from the renormalizcd Hamiltonian _ff( r ). Therefore, the 
exact ground state energy has to be within an interval 
of Ajj(oo) around the ground state energy calculated by 
the truncated S-CUT 



\E -E' \ < Ah(oo) 



(28) 



In this way, the truncation error defined by the spec- 
tral norm is no longer an abstract expression, but gives 
a practical error bound for a physical property of the 
system. 



III. ILLUSTRATIVE MODEL 

A. Double-Hard-Core-Boson 

As illustration of the formalism described above, we 
investigate the truncation error of a model of two sites 
which can be occupied by at most one particle each. To 
describe the system in second quantization, we use the 
hard-core boson language. The commutator of the asso- 
ciated annihilation and creation operators on site i and 
j is given by 



a>i,a'. 



= S.j.j ( 1 — a\a.j ] — a\a 



The Hamiltonian under study reads 

H = et + /i (ajai + aja^ ) + t (a\a2 + a\a\ 



hi 



r 21 ia\a\a2 



(l 2 



h.c. 



+ r 20 ( a\a\ + a,\a,2 ) + Va\aia 2 a2 



(29) 

(30a) 
(30b) 
(30c) 
(30d) 



Terms in the lines (|3T)c]) and (f30d) are not present in 
the bare Hamiltonian but may emerge during the flow. 
The quantity e defines the vacuum energy, fi stands for 
the chemical potential. The particle-particle interaction 
is denoted with V and t is the prefactor of the hopping 
term. The quantities T 10 , T 20 and T 21 violate the number 
of quasiparticles. They represent the non-diagonality of 
the Hamiltonian. 



B. Flow equations 



For our study, we use the particle conserving generator 
scheme rf(£) = r} pc [H'(£)], for details see Refs. 11 0, G3 
The differential equations for the coefficients of H' read 



(31a) 

(31b) 
(31c) 

(31d) 

_2p21'p20/ 

(31e) 

(31f) 
(31g) 



Since the block-band structure is conserved by ^ pc , see 
Ref. 0, HI r 20 ' stays zero during the flow unless it is 
already present in the initial Hamiltonian. As a mini- 
mal truncation scheme, we neglect the particle-particle 
interaction given by V in the following and thus the cor- 
responding contributions to d(T 21 ' and diT 20 ' . Therefore 
the only contribution to the inhomogeneity 



K (£) = (T w '(£)T 2V (l) + 8T 21 '(£)T 21 '{£)) ojoiojoa (32) 



is given by the former derivative of the particle-particle 
interaction. To make use of the possibility of calculating 
a rigorous bound for the accuracy of the ground state en- 
ergy and the maximal energy eigen-value, we choose the 
spectral norm. Since we use hard-core bosons, a\a\a\a2 
has the maximal eigen-value of unity. Equation (I20b|) 
immediately yields the truncation error 



d t e' = 


-_4P 10/ p 10 ' 


_2p20/p20/ 


dtii! = 


^plO/plO/ 


_|_2p20/p20/ 




_2p21/p21/ 


_4P 10/ p 21 ' 


dit' = 


^pl0/p21/ 


2p21/p21/ 


d e r w ' = 


-r 10 y 


-T w 't' 




_p20/pl0/ 


— 3p20'p21' 


d e r 21 ' = 


-r 21 y 


+T 2V t , 




+2r 10 't' 


_l_4P 10/ p20' 




_p21/ y / 


_ T W'y/ 


d e r 20 ' = 


-2r 20 ' M ' 


-r 20, v' 


d t V' = 


16r 10/p21/ 


_j_ g 2 1 f 2 1 / 



Ah(£) 



T 10 '(£')T 21 '{£') + 8(r 21 ') 2 {£') 



d£'. (33) 



Due to the small number of couplings, we are able to 
calculate the quantities H' and H" directly. To this end, 
we need to calculate H {£) by transforming it like an ob- 
servable under the flow of r/'(£) following Eq. We 
stress that H is transformed in this way by an exact uni- 
tary transformation without any truncations. We obtain 
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Figure 2: The vacuum energy e'(£) of H'(£) converges to the 
ground state energy E' (oo) of H'(oo). Due to truncation 
errors, the latter starts to deviate from the true ground state 
energy Eq when the flow sets in at £ = 0. The calculation is 
carried out for ^ (b) = 2,t (b) = 1 and r 10(b) = 1. 



the set of differential equations 



d e e = 


_4r 10 'r 10 — 


9p20/p20 


(34a) 


%i = 


r 10, (4r 10 


-2r 21 ) + 2r 20 T 20 






-2r 21 '(r 10 h 


~r 21 ) 


(34b) 


d f t = 


_2pl0/p21 _ 


2r 21 '(r 10 + r 21 ) 


(34c) 


d e r w = 


-r 10 '( M + * 


_|_ p20^ _ p21/p20 






_p20//pl0 _ 


-r 21 ) 


(34d) 


ftr 21 = 


r 10 '(2r 20 


+ 2t- v) + r 20/ (2r 10 + 


r 21 ') 




+r 21 '(- M 4 


■ t-v + r 20 ) 


(34e) 


ftr 20 = 


_ 2r 10/p21 + 


r 20/ (-2 y u-v A ) + 2r 21 'r 


10 (34f) 


d t V = 


8 pl0,p21 + 


8r 21 '(r 10 + r 21 ). 


(34g) 



Thereby we are able to calculate exactly as refer- 

ence to estimate the quality of the truncation error Ah 
which yields an upper bound to One should no- 

tice that the conservation of the bandstructure does not 
hold for H since the generator depends on H'. 



C. Results 

For our calculations, we used fi^ = 2 and = 1 as 
initial conditions for H<- h \ while e<^\ r 20 ( b ) and T 21 ( b ) 
are chosen to be zero. The initial non-diagonality F 10 ( b ) 
is used to control the degree of necessary transformation. 
For small values of T 10 , the Hamiltonian is close to diag- 
onality and therefore only slightly changed by the CUT. 
A large initial non-diagonality on the other hand requires 
intensive re-ordering processes in which truncation errors 
are important. 

Figure[2]shows the generic behavior of the ground state 
energy E' of H'^ and the vacuum energy e under the 



truncated flow. Truncation errors have a noticable im- 
pact due to the (large) non-diagonality T 10 . In an early 
stage of the flow (£ < 0.5), E (£) starts to depart from 
Eq and remains constant for the rest of the flow. The 
vacuum energy e' converges rapidly to the ground state 
energy E' (oo) of H'(oo). 

The truncation error Ah and the spectral distance 
1 1 #"11 both saturate in the course of the truncated flow, 
see Fig. [3k. and b, and they show a strong monotonic de- 
pendence of the non-diagonality. The spectral distance 
||i?"|| is bounded by the truncation error as it has to be. 
Furthermore, the truncation error turns out to be also a 
good approximation for the spectral norm. This can be 
seen in Fig. [3] (panel c). Their difference is insignificant 

for small non-diagonalities and even for very large ones 
( r io(b) = 10 ) it takcs Qnly 6% 

The influence of truncation on the spectrum of H' 
is studied by the difference of the ground state energy 
AEq = Eq — E' Q and by the difference of the energy of the 
highest excited level AE max = E max — E' max compared 
to the values of the initial Hamiltonian Hq. Both quanti- 
ties stay clearly below the spectral distance, but differ in 
magnitude. For small values of r 10 ( b ), AEq is negligible 
and rises only up to 3.5% of | |iT"W 1 1 for r 10 ( b ) = 10. By 
contrast, AE max is nearly identical to ||lf"'- r - l || for low 
non-diagonality and remains close to ||i/"'- r ^|| even for 
high values of r 10 ( b ). Thus the major impact of trunca- 
tion occurs for the highest excited level. 

This can be understood on recalling that the minimal 
truncation of the density-density interaction expressed 
by V means neglecting the energy correction for the dou- 
bly occupied state and to approximate it by Eq + 2/i. 
Because this is precisly the highest eigen-state for van- 
ishing non-diagonality r 10 ( b ), it is directly affected by 
the truncation. In contrast, the ground state is only in- 
fluenced indirectly by inaccuracies for the higher levels. 
Therefore, the low energy properties can be characterized 
by S-CUT very accurately despite of truncation, whereas 
larger inaccuracies occur at high energies. 

In summary, the truncation error Ah defined in Eq. 



(|20b[) is illustrated as an upper bound for the spectral dis- 
tance ||#"|| and for the errors of Eq and E max . We stress 
the good agreement of Ajj, \\H"\\ and AE mllx . The er- 
ror of ground state energy, however, is significantly lower 
than the bound given by ||J?"||. This can be explained as 
a particularity of the minimal truncation scheme which 
affects mainly high energies. We highlight that the trun- 
cation error measures the non-unitarity of the complete 
transformation acting on the whole Hilbcrt space. 



IV. EXTENDED SYSTEM 

A. Dimerized spin-^-chain 

In the previous section, we studied the truncation error 
for a zero-dimensional illustrative model. To illustrate 
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Figure 3: 

Truncation error Ah, spectral distance deviation 
A_B max of the highest and lowest energy eigen- value AEo of 
the truncated Hamiltonian for /^ b ' = 2 and = 1. Panels a 
and b: dependence of the flow parameter I for a moderate 
non-diagonality r 10 ' b ' = 1 (a) and for large non-diagonality 
pio(b) _ Panel c: dependence of the renormalized 

quantities on the initial non-diagonality r 10 ( b '. 




I 



Figure 4: Schematic representation of the dimerized spin 
chain. Dark bonds stand for the coupling J between two 
S = Y 2 spins forming a dimer, light bonds denote the vari- 
able inter-dimer coupling AJ. The dashed line indicates an 
axis of reflection symmetry. 



the applicability of the truncation error for more rele- 
vant models, we extend our analysis to a one-dimensional 
system. Furthermore, this allows us to calculate and to 
compare the truncation errors of more complex real-space 
truncation schemes. As we will see, we have to use the 
triangle inequality again to arrive at error bounds in ex- 
tended systems. 

The model studied is the one-dimensional dimerized 
antifcrromagnetic spin S = V 2 Heiscnbcrg chain with 



the Hamiltonian 

H = jJ2( S r-S? + ASf -S£ +1 ) , J>0. (35) 

r 

In this notation, stands for the operator of the left 
spin in the dimer on position r and for the right spin, 
see Fig. @] The parameter < A < 1 denotes the relative 
strength of interdimer coupling. It is used as a control 
parameter similar to T 10 in Sect. Mil 

In the limit of A = 0, the system consists of isolated 
dimers. The ground state is given by a product state of 
singlets with a ground state energy of = —jJ per 
dimer. The local 5 = 1 excitations form equidistant 
spectrum with an increment of J. 

For rising interdimer coupling, the excitations can be 
described by gapped spin S = 1 quasiparticles called 
triplons fl6l ] . They can be seen as triplets with magnetic 
polarization cloud. In the limit A = 1, the g ap closes and 
the correlations decay algebraically [3ll - l3a |. For A = 1, 
the system is the well-known homogeneous spin chain ex- 
actly solved 1931 by Bethe by what was henceforth 
called Bethe ansatz. The ground state energy for the in- 
finite chain was calculated by Hulthen [38| and takes the 
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value 1/2 — 2 In 2 per dimer and J. 

With respect to the limit of isolated dimers, we choose 
for dimer r a local basis of singlet/triplet states 



\s) r = -j= (in) - in)) 


(36a) 


\X) T = ^ (Itt) - III)) = 4, rl*) r 


(36b) 


\y) T = ^= (in) + in)) = ti r \ s ) r 


(36c) 


\*) r = ^ (in) + m» = 4» P 


(36d) 


as in the bond operator representation f39j-[4l| . We define 
the reference state as the product state of singlets on each 
site 


10) :=®|*> P . 

r 


(37) 


The triplet operators arc defined by 




4, r : = l«)r ( S lr 


(38a) 


t a ,r ■= \ s ) r Hr' 


(38b) 



Non-physical artifacts (e.g. states with two triplets on 
one dimer) are excluded. Therefore, the triplet operators 
obey the hard-core algebra 



t V 



8 rs 8 al3 1 - }j\ r t 



7,r 



S f t 



(39) 



The normal-ordered products of triplet operators (mono- 
mials) together with the identity 1 form the basis {A{\ 
for all operators on the lattice. In this notation, the 
Hamiltonian reads 

r 

+ (fa,r + ^a,r) {^a,r+l + ^a,r+l) 



By the CUT the triplet states are mapped to re- 
normalized S — 1 excitations (triplons). 

The Hamiltonian (|40D has three different symmetries 
that can be used for simplification of the calculation as 
mentioned in Sect. II PI 

(i) The Hamiltonian (|40p is self-adjoint. Although this 
is not a symmetry in the strict sense of the word, it im- 
plies an additional constraint for the coefficients of H . 
Since quasiparticle creating and annihilating terms have 



to occur in pairs in any Hamiltonian, one of them can 
be chosen as representative for the pair. Thereby the 
number of coefficients to be tracked is reduced. 

(ii) The Hamiltonian (|40| shares the reflection sym- 
metry of the chain r — > —r, see Fig. 21 In addition, 
all left-spin and right-spin operators have to be swapped 
implying t a r -f-^ — t a r in the triplon notation J45j. 

(iii) The Hamiltonian ((35)) is invariant under SU(2) ro- 
tations in spin space. Due to this invariance, the Hamil- 
tonian written in the triplet algebra (|4"0"f can be decom- 
posed into symmetric combinations of terms. The terms 
in a symmetric combination differ only by permutations 
of triplet polarizations up to a sign factor. The super- 
operator ^ G to build the symmetric combination of a 
maximally asymmetric representative reads 



E = EE = ( 1 +^ 



1 + 5 cyc + £ c 2 y c ) , (41) 



where the cyclic permutations of triplet polarizations is 
denoted by S cyc and the exchange of the polarizations x 
and y with a negative sign factor for each triplet operator 
reads S xy . 



B. Real-space truncation 

For an extended system, the omission of processes that 
create or annihilate more than TV triplons as mentioned 
in Sect. II CI can be insufficient because the number of 
remaining terms is still infinite. For example, even by 
restricting to processes of at most two triplons, an infi- 
nite number of independent terms varying by range can 
emerge in the course of the flow. Since an energy gap im- 
plies a finite correlation length, it is an adequate choice 
to neglect all processes that exceed a given range d max - 
In a one dimensional system, the range can easily be de- 
fined as the distance of the rightmost and leftmost triplet 
operator in a term. 

In particular, it has turned out to be advantageous 
to use a combination of both the quasiparticle and the 
range criterion as truncation scheme [TO, [Hj]. For the 
most important processes of low quasiparticle number, 
e.g., the hopping of triplons, a long range is allowed for to 
preserve most of the relevant physics. For more complex 
processes of more quasiparticles only a shorter range can 
be considered because their number increases much more 
steeply with range. But since more quasiparticles are 
required for such terms to become active the reduced 
range does not need to imply a reduced accuracy. 

In view of the above considerations we classify terms 
by the sum n of created and annihilated quasiparticles. 
For each value of n a specific maximal range d n is de- 
fined. The complete truncation scheme can be written 
as d = (d,2, ^3 . . . d,2N) where at most N quasiparticles 
may be created or annihilated. No maximal range needs 
to be specified for n = 1 because those terms always 
have range zero and their number is, due to translation 
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symmetry, restricted to the six local creation and anni- 
hilation operators. Furthermore, without magnetic field 
no single annihilation or creation of triplons takes place 
due to the conservation of the total spin. 



C. Triangle inequality 

To calculate the truncation error according to Eq. 
(|20b|) , the norm of the inhomogcneity 



k(£) =J2 K *MAi 



(42) 



has to be calculated. Its coefficients K,A t are obtained 
using Eq. (| 13[) from the coefficients of H' by evaluating 
the terms of the commutator that are discarded in the 
truncation scheme. The precise calculation of \\k\\ is not 
feasible for large systems because the effort to calculate 
the maximal eigen-value is too large. Since we are only 
interested in an upper bound, we apply the triangle in- 
equality again to reach 



<5>*(*)lll^ll- (43) 

Ai 



In this way, we define an upper bound Ah or the trun- 
cation error 



A H {£) 



E 



\K Ai {£')\^'\\M\> A h {£). 



(44) 



Therefore, Ah is also an upper bound for AEq and 
Ai? max , although it is less strict than Ah- 

Recall that our basis operators Ai are normal-orderd 
products of hard-core-boson creation and annihilation 
operators. It turns out that in this case the spectral norm 
1 1 Aj| | can be calculated easily since the product Ai is 
already diagonal with respect to the dimer eigcn-states. 
It is a product of local triplon density terms. Thus A4' A4 
has the eigen-value unity for all configurations having 
triplons with the polarization of the local trip lon density 
operators on each site of the cluster of Ai |46( ■ For differ- 
ent occupations, A^ Ai has the eigen-value zero. There- 
fore the spectral norm for all elements of the operator 
basis {Ai} is unity. 



D. Optimization using symmetries 



As pointed out in Sect. IIP! exploiting symmetries re- 
duces the computational effort for the S-CUT method 
significantly. In addition, symmetries can be used to re- 
duce the bound A. To see this, we write A using the basis 



of representatives {Ci} of operator monomials 



K c z 



df 



(45a) 



<E / di> - 



G 



=:A H (£). (45b) 



We use the triangle inequality again to decompose 
I Eg ^11 m t° the weight wet which denote the number 
of terms generated by multiplied by the norm of the 
representative. This yields 



A„{1) 



\K C ^')\d£'\wcA\Gi\\. 



(46) 



If, however, we avoid the triangle inequality for the 
fully symmetric combination of basis operators we obtain 
a better, stricter upper bound. This enters Eq. (|4"6]) by 
replacing wet by a reduced effective weight Wq. . In order 
to determine this reduction, we have to calculate the re- 
lation between the norm of the fully symmetric combina- 
tion and the norm of a single representative analytically 



E^ 



\Ci 



(47) 



The effective weight factor Wq. allows us to modify Eq. 

(|4"6")l to calculate the improved truncation bound A G . In 
presence of multiple symmetries, it is possible to improve 
the weight with respect to only some selected symmetries. 
The other symmetries can be treated by using the trian- 
gle inequality as done in Eq. (|46[) . They continue to enter 
the combined weight by an additional factor. 

As an example, we consider the self-adjointness to 
reach the improved truncation bound A^. To determine 
its weight, it is useful to decompose the action of a rep- 
resentative 



C=\c l )(c a \ Q§ 1 

r^cluster 



(48) 



into the non-trivial action on its cluster and into the 
identity on the rest of the system. Since C is a normal- 
ordered product of hard-core-boson operators, the action 
on its cluster is given by only one non-vanishing matrix 
element. 

To build the corresponding symmetric combination 
C , we have to distinguish two cases: 

(i) |ci) = \c2) If C is self-adjoint, it is already a 
symmetric combination and thushas a weight factor of 
unity in both cases, i.e., using and not using the self- 
adjointness. 

(ii) \c\) 7^ \c2) In this case, the symmetric combina- 
tion for C reads 

^C=(l + l)c=(| Cl )( C2 | + |c 2 )( Cl |) (g) 1. 

G r^cluster 

(49) 
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Due to our choice of basis, both |ci) and |c 2 ) are eigen- 
states of local triplet density operators and therefore or- 
thogonal, which implies the maximal eigen-value of 1 for 
|ci) (c 2 \ + \c 2 ) (ci|. Hence || \c x ) (c 2 | || = 1 and || | Cl ) (c 2 | + 
\ c 2) (ci \ || = 1 which implies a gain of a factor of 2 if the 
triangle inequality is avoided which would have led to 
II |ci)(c 2 | + |c 2 )(c 1 |||<|||c 1 )(c 2 ||| + |||c 2 )(ci|||=2. As 
an example, we consider the representative t x 1 t x \i y2 - 
On its cluster {1, 2} it stands for a transition from |xis 2 ) 
to \xiy 2 ). The matrix representation for the action of its 
symmetric combination 

^U^U = (\^V2) (X1S 2 \ + \X1S 2 ) (x m \) (g) 1 
G r£{l, 2 } 

(50) 

has zero matrix elements except for a 2 x 2 block with 
the eigen- values -1 and +1. 

In conclusion, using sclf-adjointness to gain an effective 
weight saves a factor of 2 for non-symmetric terms. 

This effective weight can be improved further by con- 
sidering spin symmetry. The complete symmetry group 
with respect to all permutations of triplet polarizations 
can be decomposed into the subgroup of cyclic permu- 
tations and the subgroup of the transposition of x and 
y triplets S xy plus the identity. For simplicity, we con- 
centrate on the latter one only and calculate the trunca- 
tion bound A H xy . The matrix associated with the J^g ^ 
action on its cluster can have up to four non-vanishing 
matrix elements , i.e., four states of the cluster have to 
be taken into account. The representatives can be clas- 
sified by the specific action of ^ G needed to obtain the 
corresponding fully symmetric combination: 

= 1- These highly symmetric representatives (e.g. 
t z 1 t z jj are invariant under both A and S xy . They have 
weight w*' xy unity. 

Y2g = 1 + A: These representatives are not self- 
adjoint, but invariant under either transposition of x and 
y (e.g. t z jj or under the combination of transposing and 
adjunction (e.g. t x ^t y i ). As discussed previously, the 
weight w'> xy takes the value of 1 instead of 2. 

^2 G = 1 + S xy : For representatives that are self- 
adjoint, but not invariant under transposition (e.g. 
t x 1 t x l ), the weight w^ xy is reduced to the value 1 in- 
stead of 2 as well. 

J2g = (^^^) (■"■ ^ x v) ; ^ n ^ ms asymmetric case, 
the norm of the symmetric combination can be either one 

( e -g- 4,i*»,a) or (e.g. 4,i)- Therefore, u>t.*» = ^2 
can be used instead of w = 4 . 

The weights can be improved by exploiting further 
symmetries, e.g., cyclic spin permutations or reflection 
symmetry. But the complexity of the necessary case- 
by-case analysis rises considerably. Especially the use of 
point group symmetries of the lattice is complicated be- 
cause the cluster of monomials linked by a point group 
do not need to be identical. Thus the inclusion of point 
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Figure 5: Renormalized truncation bound per dimer Ah and 
reduced bounds (exploiting self-adjointness) and A H xy 
(exploiting self-adjointness and xy symmetry) vs. the inter- 
dimer coupling A for the ground state generator using the 
truncation scheme d = (8, 6, 6, 3, 3). 
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Figure 6: Renormalized reduced truncation bounds per dimer 
A^ xy vs. interdimer coupling A for ground state and particle 
conserving generator using various truncation schemes. 

group symmetries in the calculation of the bounds is be- 
yond the scope of this article. 



E. Results for an infinite system 

Figure [S] shows the bound Ah defined in Eq. ([4T)f for 
the infinite dimerized Heisenberg chain using the ground 
state generator scheme 7/ gs [l9| defined in (0 and the 
truncation scheme d = (8,6,6,3,3) as function of the 
interdimer coupling A. For rising values of A, the trunca- 
tion bound increases drastically similar to the behavior 
observed for the illustrative model. But it attains a sig- 
nificantly higher absolute value finally. 

For weak interdimer coupling A < 0.1, the error bound 
for the ground state energy given by A is useful as a 
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rj d | Ag 1 A H A h A H xy 

77p C (6,3,3) 0.01739 439.37 219.93 156.74 

?7 PC (8,6,6,3,3) 0.00032 49926.1 24974.8 17660.2 

%b (6,3,3) 0.02675 23.51 11.75 8.41 

?7g S (8,6,6,3,3) 0.00915 149.31 74.68 52.98 

Vzs (9,7,7,4,4) 0.00948 169.32 84.67 60.04 

Table I: Numerical values for the renormalized truncation 
bound per dimer Ah and reduced bounds A H , A H xy for A = 1 
for the generators and the truncation schemes used in Fig. [5] 
The inaccuracies of the ground state energy |A_Eq| are calcu- 
lated with respect to the analytical result [38q . 

rigorous bound. But this estimate becomes inappropriate 
for medium and strong coupling since the A grows rapidly 
to the same magnitude as E and beyond so that it does 
no longer represent a meaningful bound. 

For A = 1 the analytical result can be used as refer- 
ence to determine the error of the ground state energy 
per dimer |Ai?o| determined as the renormalized vacuum 
energy of the CUT, see Tab. [I] This relative error is only 
1.04%. In contrast, the truncation bound Ah exceeds 
the ground state energy by several orders of magnitude. 

This discrepancy stems from the fact that the trunca- 
tion error Ah measures truncation effects of the entire 
transformation, not only of the inaccuracies of ground 
state energy in particular. The extensive use of the tri- 
angle inequality to calculate the truncation bound A en- 
hances this difference additionally, although it can be 
reduced exploiting symmetry. The interesting question 
which effect dominates is postponed to the study of a 
finite system below where all quantities are numerically 
accessible. 

The use of the adjunction symmetry reduces the bound 
by a factor of about two, and the use of xy symmetry re- 
duces it by an additional factor of about \/2. In both 
cases, exploiting the symmetry pays in decreasing the 
bound close to the optimum which can be achieved for 
terms with the lowest symmetry. Since terms with low 
symmetry occur much more frequently than symmetric 
ones, the error bound is reduced efficiently by exploiting 
symmetries. However, the gain achieved in this way can 
not overcome the tremendous factor (orders of magni- 
tude) between the upper bound and to the error of the 
ground state energy. 

In the extended system, various truncation schemes of 
different quality and computational effort can be used. 
Figure [5] shows the truncation bound A H xy for different 
truncations and generator schemes. Besides the parti- 
cle conserving generator rj pc the ground state decoupling 
generator ^ gs is applied, for details see Ref. [r| The nu- 
merical values for A = 1 including the difference to the 
exact ground state energy are given in Tab. [I] 

In general, the particle conserving generator scheme 
Vpc yields significantly higher truncation bounds than 
77g S , although the inaccuracies of the calculated ground 
state energies are lower (using equal truncation schemes 
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Figure 7: (Color online) Renormalized truncation bound per 
dimer A H xy , exact renormalized truncation error per dimer 
Ah and error of ground state energy per dimer |Ai?o| vs. in- 
terdimer coupling A for the ground state generator considering 
at most N triplon operators per term. 



d). This is due to the fact that rj pc performs a more 
comprehensive reordering of the quasiparticlc subspaccs 
and includes much more terms than r} gs . This implies 
a higher impact of truncation errors. Much more terms 
emerge in the evaluation of the commutator and have 
to be incorporated in n resulting in a larger differential 
equation system with much more coefficients. 

The dependence on the truncation scheme is more com- 
plex. For weak coupling, looser truncation schemes im- 
ply lower truncation error bounds than stricter schemes. 
This is what one expects naively since the inclusion of 
more and more terms, hence a less strict truncation, 
should describe the system better and better. Thus the 
truncation error should decrease. 

But this relation can be inverted for strong coupling 
although the looser schemes reproduce the analytical re- 
sult for the ground state energy with much higher ac- 
curacy for both generator schemes as seen for (6, 3, 3) 
and (8,6,6,3,3). We call this phenomenon the trunca- 
tion paradoxon because it seems to be paradoxical at first 
glance. We stress that it does not represent a logical con- 
tradiction but only a counterintuitive behaviour. 

On second thought, one realizes that the use of a looser 
truncation scheme increases the number of terms in n 
drastically. For instance, an increase from 30.972 for 
(6,3,3) to 16.777.215 representatives for (8,6,6,3,3) is 
found in the fj pc scheme. Irrespective of the consequences 
to the spectral norm of k, this massive increase of terms 
has a big impact on the bound A which relies on the tri- 
angle inequality to bound each of these terms separately. 
Hence the bound is so large simply because it is very 
loose. 
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N 


\AE \ 


A H A h 


aL 


At. 

IV H 


2 


0.04109 


0.1556 20.34 


10.17 


7.28 


3 


0.01709 


0.0964 41.87 


20.93 


14.86 


4 


0.00465 


0.0768 26.78 


13.39 


9.48 



Table II: Error of ground state energy per dimer |A_Eo|, exact 
renormalized truncation error per dimer Ah and renormalized 
truncation bounds per dimer A H , A H and A H xy for A = 1 
using the ground state generator scheme rj ss considering at 
most N triplon operators per term. 



F. Results for a finite extended system 

The question arises whether the truncation paradox is 
caused by the extensive use of triangle inequality to de- 
termine the bound A, or whether it is an intrinsic char- 
acteristics of the truncation error A itself. To investigate 
this issue, it is necessary to determine the exact trunca- 
tion error A without use of the triangle inequality. Thus 
an exact diagonalization of K is required. This restricts 
us to the investigation of a finite chain segment. 

In the following, we study the periodic, dimcrized 
Heisenbcrg chain consisting of five dimcrs. In view of 
the small size of the system, we use the maximal number 
N of interacting triplons as only truncation criterion. No 
extensions are considered. 

Figure [7] shows the truncation bound A H xy , the exact 
truncation error A and the deviation of the ground state 
energy |Ai?o| vs. A. The values for A = 1 are given in 
Tab. |nj It turns out that both the calculation of the exact 
truncation error Ah and the extensive use of the triangle 
inequality contribute to the very large values of A H xy . 
For small values of A, the large value of Ah dominates 
the truncation bound A H xy while for large values of A, the 
approximation using the triangle inequality to bound the 
very many arising terms contributes most. 

In a direct comparision, the exact truncation error Ah 
is overestimated by the bound A H xy using the triangle in- 
equality by two orders of magnitudes, even though sym- 
metries are used. Nevertheless, Ah is still considerably 
higher than the deviation of ground state energy. In par- 
ticular, it exceeds |Ai?o| by several orders of magnitudes 
for small A. Here we have to keep in mind that the trun- 
cation error does not only measure the inaccuracies of 
ground state energy, but the effect of truncation to the 
entire transformation of the Hamiltonian. 

In contrast to the double hard-core boson model, no 
deviations in the maximal energy eigen-value were ob- 
served (not shown). This is explained by the fact that 
the fully polarized state is still an exact eigen-state of the 
Hamiltonian (|35[) . 

For high values of A, the truncation paradox occurs 
again because the truncation error bound A H xy for the 
three-triplon and the four-triplon truncation exceeds the 
truncation error bound of the two-triplon truncation. 
The exact truncation error Ah does not display any para- 



doxical behaviour. Hence we conclude that the extensive 
use of the triangle inequality is at the basis of the trun- 
cation paradox. 

V. SUMMARY 

In this work, we presented a mathematically rigorous 
framework to bound effects of truncation in self-similar 
continuous unitary transformations a priori. The dif- 
ference H" = H — H' between a unitarily transformed 
Hamiltonian H and the Hamiltonian H' obtained from 
the truncated calculation is captured by an inhomoge- 
neous flow equation depending only on the truncated 
terms. We defined the scalar truncation error A by the 
norm of the truncated terms ||k||. It provides an upper 
bound for the norm of the difference 1 1 H" 1 1 . A completely 
analogous bound is derived for observables as well. 

Using the spectral norm, the truncation error implies 
an upper bound for the deviation of the minimal and 
maximal eigen-value of the observable under study, which 
is caused by truncation. In particular, we derived a rig- 
orous a priori bound for the error of ground state energy. 

The analysis of the double hard-core boson model 
showed that the norm of the difference H" is bounded 
and approximated very well by the truncation error. De- 
spite the large difference to |A-Eo|j the truncation error 
Ah provided a good measure for the deviation in the 
highest excited level |A.E max |. This could be understood 
by the special feature of the truncation scheme that pri- 
marily affected the highest excited level. 

For practical use in extended systems, an upper bound 
A for the truncation error is calculated using the triangle 
inequality. The direct calculation of ||k|| is not feasible 
- even impossible for infinite systems - because it would 
require an exact diagonalization in the Hilbcrt space of 
the entire system. 

In both systems studied, the double hard-core boson 
model and the extended dimcrized spin chain , the bound 
provided by truncation error Ah turned out to overesti- 
mate the actual inaccuracies of the ground state energy 
significantly. This is an inevitable consequence of the 
fact that the truncation error is a measure for the non- 
unitarity of the whole transformation. Therefore it is sen- 
sitive to distortions of all eigen-values and eigen-states. 

A comparison of various bounds for finite dimcrized 
spin chain segments showed that the truncation bound 
Ah overestimates the real truncation error Ah by orders 
of magnitudes. Furthermore, a truncation paradoxon was 
observed: Looser, i.e., better, truncation schemes implied 
higher truncation bounds. In the finite system studied 
using the ground state generator, this paradoxon did not 
occur for the exact truncation errors. 

The truncation bound can be improved efficiently ex- 
ploiting symmetries of the Hamiltonian and its her- 
mitecity. By using the transposition of the triplon po- 
larizations x and y and the hermitecity, we were able to 
reduce the bound by a factor w 2y/2. To overcome the 
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L 


A = 0.5 


A = 1.0 


3 


-0.7821366540 


-0.9342579563 


4 


-0.7778813863 


-0.9019258071 


5 


-0.7770436419 


-0.8836215559 


6 


-0.7768700368 


-0.8745235739 


7 


-0.7768649893 


-0.8740959679 


8 


-0.7768649893 


-0.8740959679 


9 


-0.7768649893 


-0.8740959679 


10 


-0.7768649893 


-0.8740959679 



Table III: Ground state energies per dimer Eo for the finite 
dimerized Heisenberg chain obtained by S-CUT using a real- 
space truncation scheme for systems of different numbers of 
dimers n and different interdimer couplings A. The maximum 
untruncated range is set to d max = 2, the corresponding ef- 
fective system size La n is 7. Beyond this size, the calculated 
quantities become independent of the system size. 

high difference to the truncation error A, much more so- 
phisticated approximations would be needed. Additional 
symmetries available are the cyclic spin permutation, the 
reflection symmetry, and the translation symmetry. We 
do not see a way to exploit the powerful translation sym- 
metry completely for the improvement of the bounds be- 
cause this requires finding bounds for operators in very 
large or infinite systems. But larger subsets of terms of 
restricted range could indeed be analyzed on larger clus- 
ters. The calculation of H" in an extended system is left 
for further investigation. 

The analysis presented here provides a better under- 
standing of truncation errors in continuous unitary trans- 
formations. We are confident that it serves as seed for 
stricter a priori error bounds in the future. 
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Appendix A: Effective system size 

In Sect. Ill Al we mentioned that using a real-space trun- 
cation scheme implies an effective system size £fi n that 
disguises the differences between infinite and finite peri- 
odic systems with a size of at least L^ n . The coefficients 
of representatives calculated by a truncated S-CUT be- 
come independent of the physical system size. At first 
glance, this seems to be puzzling. But one has to look in 
detail where the finite system size enters in the calcula- 
tion of the flow equation. 




Figure 8: (Color online) Schematic representation of the gen- 
eration of a term which complies with the truncation scheme 
with d max = 2 only due to the finiteness of the systeme with 
size I = 6. Following the hard-core-algebra, the commuta- 
tor of the terms t\ x t\ x and t3 tX t\ x (black) generates several 
terms, see Eq. (I Al|) . with range 4 that are truncated (or- 
ange/gray and dashed), but the term ti tX t\ x is not truncated 
due to a ,wrap-around' effect (blue/gray). 



For clarity, we focus on a one-dimensional system, al- 
though the argument can be generalized to higher dimen- 
sions. Let d max be the maximal range of the untruncated 
terms and I the systems size. In order to keep the trans- 
lational invariance of the infinite system also for a finite 
chain, periodic boundary conditions must be chosen. By 
evaluation of the flow equation ([2]), the range of the terms 
stemming from the commutator [>](£), H(i)] becomes im- 
portant. 

For an infinite system, the issue is clear because the 
maximum range of the commutator of two processes of 
range d\ and di is given by the sum d comm = d\ + d-2 
d c = di + d%. In a periodic system the range of a term is 
a more complex issue. The cluster of a term ambiguous 
because each site may be chosen to be the leftmost site.. 
As an example, the range of the term t\ x t\ x in a system 
of six dimers can be defined to be either 4 or 2, see Fig. [8] 
To remove this ambiguity the most plausible prescription 
is to define the range of a term to be the minimum range 
of all possible choice of the leftmost site. 

If the maximal range in the generator and in the Hamil- 
tonian is rf max the maximal range of terms stemming from 
the commutator is given by d C omm,max = 2d max . Due to 
the locality of the algebra, the clusters of both terms 
have to share at least one common site as can be seen in 
Fig. [5] So for smaller systems L < d CO mm,max the result 
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of the commutation is strongly influenced by the finite 
size. If the system is larger than d C omm,max, the terms 
from the commutator are identical to those in the infi- 
nite system. But even then a wrap-around can influence 
the range which is actually attributed to a term result- 
ing from the commutator. As an example we consider 
the commutator 



b l,x L 3,xi b 3,x t/ 5,x 



b l,x b 5,x 



l.x 3.x 3.x 5,x 



^ b l,y b 3,y b 3,y b 5,y + b l,z b 3,z b 3,z b 5,z- 

In the infinite system, each of the terms on the right- 
hand side of (|Aip would be truncated since their range 
4 exceeds the maximal range set to 2. But in the peri- 
odic system, the range of the contribution — 1\ x t\ x can 
be lower due to a ,wrap-around' if the system is small 
enough, see Fig. [5J 

In general, this anomalous range assignment due to a 
wrap-around can happen if the relation 



L < d, 



comm. max 



^max — 3flxi 



is fullfillcd. As a result one obtains the same differential 
equations for the representatives for all systems with a 
size of at least 



I/fin = 3d n 



(A3) 



(Al) 

This included the infinite system. Hence Lfi n defines the 
effective size of a system treated with the maximal trun- 
cation range d max . 

For a numerical illustration, we examine the depen- 
dency of the ground state energy per dimer for the dimer- 
ized Heisenberg chain introduced in Sect.[lV] We use the 
real-space truncation scheme with a maximal truncation 
length d max = 2 and no restrictions of the number of in- 
teracting triplons. This implies an effective system size 
Lg n of 7 dimers according to (|A3[) ■ The results are given 
in Tab. [21 Clearly, the results are numerically identical 

(A2) 

for larger systems as predicted. 
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